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ABSTRACT 

We investigate the distribution of massive black holes (MBHs) in the Virgo cluster. 
Observations suggest that AGN activity is widespread in massive galaxies (A/* > 
1O 1O M0), while at lower galaxy masses star clusters are more abundant, which might 
imply a limited presence of central black holes in these galaxy-mass regimes. We ex- 
plore if this possible threshold in MBH hosting, is linked to nature, nurture, or a 
mixture of both. The nature scenario arises naturally in hierarchical cosmologies, as 
MBH formation mechanisms typically are efficient in biased systems, which would later 
evolve into massive galaxies. Nurture, in the guise of MBH ejections following MBH 
mergers, provides an additional mechanism that is more effective for low mass, satel- 
lite galaxies. The combination of inefficient formation, and lower retention of MBHs, 
leads to the natural explanation of the distribution of compact massive objects in Virgo 
galaxies. If MBHs arrive to the correlation with the host mass and velocity disper- 
sion during merger-triggered accretion episodes, sustained tidal stripping of the host 
galaxies creates a population of MBHs which lie above the expected scaling between 
the holes and their host mass, suggesting a possible environmental dependence. 



1 INTRODUCTION 



The nearby Virgo cluster is a perfect laboratory to investi- 
gate the evolution of galaxies in a dense environment. Re- 
cently, observations of a large sample of galaxies suggested 
that the properties of nuclei, either quiescent or active, in 
Virgo galaxies are strongly mass-dependent. This latter find- 
ing is in very good agreement with the general trend found 
also in the SDSS, where very few A GN are found in galax- 
ies with stellar mas s M* < 10 10 M Q (|Kauffmann et aLlBooll; 
iKewlev et alll2006h . 

Decarli et al. (2007) analyzed nuclear activity in late 
type galaxies in the Virgo cluster. They conclude, quite 
remarkably, that at galaxy mas (0 M gal > 10 10 5 M the 
AGN fraction is unity. As a central black hole is a neces- 
sary condition for AGN activity, we conclude that the black 
hole occupation fraction (BHOF) must be unity as well. 
ICote et all l|2006l ). Wehner & Harris (2006) and Ferrarese 
et al. (2006) find that, below M, ~ 10 10 M , Virgo galaxies 
exhibit nuclear star clusters, whose mass scales with M* in 
the same fashion as those of the massive black holes detected 
in br i ghter galaxies (IMagorrian et al ll99Sl ; lMarconi fc Hunt] 
120031 : iHaring fc RudboOfl T Although the existence of a nu- 
clear star cluster does not rule out a small, hidden MBH, 



1 Decarli et al. (2007) measure the dynamical mass of the galaxy 
within the optical radius, determined at the 25 th mag arcsec - 2 
isophote in the B band. This is an upper limit to the stellar mass 
of galaxies, but a lower limit to the total baryonic mass. 



it is suggestive that Ferrarese et al. (2006) conclude that 
"bright galaxies often, and perhaps always, contain super- 
massive black holes but not stellar nuclei. As one moves to 
fainter galaxies, nuclei become the dominant feature while 
MBHs might become less common and perhaps disappear 
entirely at the faint end." 

There are three interlaced sides of the intriguing story 
which appears to link stellar nuclei and MBHs: 1) under- 
stand if (and why) MBHs populate preferentially bright 
galaxies; 2) understand why stellar nuclei populate prefer- 
entially faint galaxies and 3) understand why the ratio of 
nuclear to galaxy mass is identical to the ratio of MBH to 
galaxy mass. In this paper we address the first issue. A possi- 
ble hint to explain the predominance of star clusters in small 
galaxies may come from comparing the dynamical time scale 
and the fragmentation time scale of the infalling gas. De- 
tailed calculations are necessary to test this hypothesis, but 
it can be argued that in shallow potentials gas could frag- 
ment before reaching the center of the galaxy. This is even 
more suggestive in the case of merger-induced gaseou s infall . 
Regarding the third issue, lEmsellem fc van de Venl (|2007I) 
for instance show that if galaxies have Sersic profiles, radial 
compressive forces trigger the collapse of gas in the central 
regions, and the mass of the nuclear cluster that forms is 
about 0.1% of the mass of the host galaxy. So, nuclear clus- 
ter formation and MBH feedback might produce the same 
scaling relation with galaxy mass, but it is unclear whether 
this is a coincidence or the result of a single, unexplored 
process. 
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Volonteri et al. (2007) investigate the overall distribu- 
tion of MBHs in galaxies, as a function of the host velocity 
dispersion, <r* , as traced by the host halo (cfr. Ferrarese 
2002), and find that the efficiency of MBH formation 
decreases with halo mass, and isolated dwarf galaxies 
are most likely to be devoid of a central MBH. This is a 
common feature of MBH formation models which invoke 
gas-dynamical proce sses in the high redshift Universe, either 
via d ir ect collapse ( Hachncl t fc Reesl 1 19931: lLoeb fc Rasid 



19941: Eiscnstein fc Loebl 



1995 



Bromm fc Loebl l2003f : 



Begelman et al. 20061 : Lodato fc Nataraianl 20061 . e.g.,), 
or vi a an intermediate stage of metal free Population HI 
star dCarr et al.lll984l : IMadau fc Reesll200ll : IVolonteri et~ai] 
120031 . e.g.,). The common feature is the need for deep 
potential wells at cosmic epochs when the protogalaxy 
population was dominated by minihalos (Mh < 1O 7 M0). 

A second physical phenomenon strengthens the like- 
lihood that dwarf galaxies, even if seeded with a MBH at 
early times, are now lacking one: the gravitational recoil. 
That is, the general-relativistic effect which imparts velocity 
to the center of mass of a merging MBH binary, when the 
MBHs plunge into each other's last stable orbit by emission 
of gravitational radiation. Gravitational waves carry, in 
general, a non-zero net linear momentum, which establishes 
a preferential direction for the propagation of the waves. As 
a consequence, the center of mass of the bina ry recoils in the 
opposite direction l|Redmount fc R ocs 1989), possibly caus- 
ing the ejection of MBHs from the po t ential wells of their 
host galaxies (e.g., IMadau et al. 2004 1_ Madau fc Quataert] 
2004 iMerritt et al.1 12004 IVolonteri fc Reesl 120061 : iHaimanl 



2004 : ISchnittman fc Buonannol 120071 1." The progress of 



numerical relativity is now leading to a convergence in the 
estimates of the recoil. Schwarzsch ild, i.e., non-spinning, 
black holes (e.g., iBaker et al.ll2006T ) are expected to recoil 
with velocities below 200 kms -1 , and a similar range is 
expected for black holes with low spins, or with spins 
(anti-)aligned with the orbital axis. However, when the 
spin vectors have opposite directions and are in the orbital 
plane, th e recoil velocity ca n be as la rge as a few thousands 
km s~ 1 (ICampanelli et al J l2007al : iGonzalez et all 120071 : 
ICampanelli et al.ll2007bT v 

Schnittman (2007) shows that the hierarchical nature 
of galaxy assembly implies that ejections do not lead to 
void nuclei, even in the high-recoil limit. This is weaved into 
the very nature of galaxy assembly, via a series of mergers, 
in a pattern typically dubbed "merger tree" . The botanical 
essence of the tree implies that many branches converge into 
a central trunk, so that in every generation of the tree, the 
number of galaxies decreases, while the fraction of MBHs can 
increase, even if ejections operate. This is especially true for 
"central galaxies" , which represents the trunk of the merger 
tree. 

The situation is different for satellites in a galaxy clus- 
ter: they correspond to loose branches that do not merge 
into the main trunk. Satellites are indeed galaxies that en- 
ter the dark matter halo of the cluster, but do not merge 
with the central galaxy. 

In this paper we develop simple models that describe 
the dynamical evolution of satellite galaxies in a cluster, 
focusing on the fate of their nuclei. We estimate the influence 
of MBH formation mechanisms, and of MBH ejections for 
shaping the BHOF in the Virgo cluster today. 



2 MBH FORMATION AND DYNAMICAL 
EVOLUTION 

We follow the evolution of the MBH population in a ACDM 
Universe along the history of a cluster-size halo. Our tech- 
nique and cosmological framewor k is similar to the one de- 
scribed in IVolonteri et al. I (|2003h . However we now focus 
on the merger history of halos whi ch have a total mass 
M h = 1O 15 M at z = (see also lYoo et ail I2007T ). We 
track the dynamical evolution of MBHs ab-initio and fol- 
low their assembly down to z = 0. We adopt very simplified 
assumptions on MBH formation and mass growth, based on 
a phenomenological approach. Our goal is to highlight the 
influence of two main processes, efficiency of MBH formation 
and dynamics, rather than develop an omni-comprehensive 
evolutionary model plagued by free parameters. 

Several theoretical arguments indicate that MBH for- 
mation proceed s at very high redshift, and probably in bi- 



ased halos (e.g., Madau fc Reesll200ll : lHaiman 2004; S hapiro! 



120051 : IVolonter i & Rces 2005, 120061 ) . Additional paramefc 
such as the angular momentum of the gas, its ability to cool 
and its metal enrichment, are likely to set the exact efficiency 
of MBH formation and the redshift range w hen the mecha- 
nism operates (for a thorough discussion see I Volonteri et all 
l2007f) . 

We explore here 4 simplified models: a very high ef- 
ficiency model (model I, "higheff"), where all halos with 
a formation redshift z > 5 host a MBH, an intermedi- 
ate efficiency model (model II, "mideff"), where all halos 
with a formation redshift z > 12 host a MBH, and finally 
two biased, low efficiency model, where MBHs form only 
in halos which represent density peaks with v c = 50( model 
III, "peak3") or v c — 3.5 (model IV, "peak3.5") and have 
a formation redshift z > 12. Th e choice v c = 3.5 i n 
model IV is based on the results in IVolonteri et~aH i|2003h . 
and is meant as a hard threshold to illustrate the ef- 
fect of biased formation. Our models link the formation 
of a MBH seed to the halo mass and formation redshift 
of a given galaxy only, and we do not assume here that 
MBHs ar e formed in bulges. Several models of MBH forma- 
tion (e.g..lEisenstein fc Loeblll995l:lKoushiappas et~aL 1 12004 
iBegel man et alj|2006l : lLodato fc Nataraianl I2006T I addit ion- 
ally require that MBH seed hosts have low angular momen- 
tum, thus hindering the existence of MBHs in pure disc 
galaxies. 

Our four "models" of MBH formation define how often 
and when halos are populated with MBHs. They provide the 
initial occupation fraction, and they set the minimum mass 
of MBH hosts today. The higher the redshift when MBH 
formation ceases, the larger the typical MBH host today is. 
This is a dir ect consequence of hie rarchical cold dark matter 
models fcfr. IMadau fc Reesl l200lh . 

Additionally we have to follow the dynamical evolution 
of halos and embedded MBHs all the way to 2 = in or- 
der to determine the occupation fraction of MBHs and their 



2 v a = S c (z)/a(M, z) is the number of standard deviations which 
the critical collapse overdensity represents on mass scale M. 
a{Mf l , z) is the root mean square fluctuation of the linear density 
field at redshift z, and 5 c (z) is the threshold density for collapse 
of a homogeneous spherical perturbation at redshift z. See, e.g., 
Peebles 1993 
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properties. We focus here on the effect that MBH ejections, 
namely due to gravitational waves recoils, have on the prop- 
erties of the MBH population at z = 0. The magnitude of 
the recoil depends on the mass ratio of the merging MBHs, 
the spins of the MBHs, the orbital parameters of the binary. 

First, to evaluate mass ratios, we have to model the 
mass-growth of MBHs. We base our modeling on plausible 
assumptions, supp orted by both simula tions of AGN trigger- 
ing and feedback l|Springel et al.ll2005l ). and analysis of the 
relationship betwe en MBH masses (Mbh) and the proper- 
ties of their hosts dMcLure fc Dunlodl2004l ; IWvithe fc Loebl 
120051 ). IWvithe fc Loebl (|2005h show that if the relationship 
between the mass of a MBH and t he velocity dispersion 
of the host found for loc a l galaxies ( Tremaine et al.l |2002| ; 
iFerrarese fc Merrittl I2OO0I ; iGebhardt et all 12000^ does not 
evolve with redshift, then the correlation between the masses 
of MBHs and their hosts evolves with redshift in a way 
compa tible with obse rvational results bv[McLure fc Dunloa 
HOOJ). Additionally. ISpringel et~aH i|2005F suggest that the 
ditto Mbh — o* relation is established during galaxy mergers 
that also fuel MBH accretion and form bulges. We therefore 
assume that after every merger between two galaxies with 
a mass ratio larger than 1 : 10, their MBHs attain the mass 
predicted by the Mbh — c* for each of the merging galaxies. 

Hence, although in our models the presence of a MBH 
is not uniquely coupled with bulge formation, the mass of 
a black hole is set during the same event that forms the 
host bulge. Accordingly, when a binary of MBHs merge, 
their mass ratio scales with the Mbh — 0* relation appro- 
priate for the velocity dispersions of the progenitor halos of 
the MBHs. Note that if the Mbh — o* relation scales with 
redshi ft as suggested by, e.g. IWoo et all (|2006l ); lTreu et all 
the mass ratio of merging binaries would be un- 
changed, so the occupation fraction of black holes would not 
be affected. We further assume that MBHs merge within 
the merger timescale of their host halos, which is a likely 
assumpti on for MBH binar i es formed after gas rich galaxy 
mergers (jEscala et al.ll2004l ; iDotti et al.ll200d 120071 ). 

To determine the efficiency of MBH ejections due to 
gravitational recoils, we need also information on the magni- 
tude and orientation of MBH spins at the time of the merger. 
We will express MBH spins as a function of the dimension- 
less parameter a = Jh/Jmax = c 4/GMmbh, where J h is 
the angular momentum of the black hole. 

Non-spinning MBHs, or binaries where MBH spins 
are aligned with the orbital angular momentum are ex- 
pected to recoil with velocities below 200 kms -1 . The re- 
coil velocity is largest for MBHs with large spins, when 
the spin vec tors have opposite direct i ons and are in the or- 
bital plane ([C amp anelli et al.l l2007al : iGonzalez et~ail 120071; 
ICampanelli et al.ll2007bT r Assuming that BHs at the time of 
the m erger always hav e antialigned spins in the orbital plane 
(as in IVolonterill2007l ) would provide a strict upper limit to 
the effect of the recoil. However, the configuration yielding 
the highest reco il velocities is probably ra t her uncommon, as 
pointe d out by IBogdanovic et all (|2007r i . IBogdanovic et all 
(2007) suggest that when the MBH merger happens in a gas 
rich environment, and is accompanied by accretion, the most 
likely configuration has spins aligned (or anti-aligned) with 
the orbital angular momentum, thus avoiding the highest 
recoil velocity. Conversely, in gas poor mergers there is no 
preferential spin alignment, so all spin/orbital parameters 



configurations are equally probable. We will assume in the 
following that orbital parameters and spin configuration are 
isotropically distributed, likely providing a soft upper limit 
to the strength of the recoil. 

The distribution and most probable value of MBH 
spins is observationally largely unconstrained. There 
is evidence that MBHs in some local AGN ga laxies 



do spin (Istreblvanska et alj 120051 ; IComastri et al. 120061 ; 
iBrenneman fc R eynolds 2006 |), based on ir on line profiles 



(|Milleijl2007l ; iFabian etailll989l ; [Laorlll99ll ). High spins in 
bright quasars are also indicated by the high radiative ef- 
ficiency of quasars, as ded uced from o b servations, by ap- 
plying Soltan's argument l|Soltanl 1 19821 ; IWang et all 120061 . 
and references therein ). However, neither observation firmly 
establishes that most MBHs have large spins, although 
there are theoretical argume nts to expect so ( Mo derski et al.l 
1 19981 ; IVolonteri et al.l 120051) as spin-up is a natural conse- 
quence of prolonged disc-mode accretion for any hole that 
has (for instance) doubled its mass by c apturing mate- 
rial with con st ant angular mome ntum axis (|Bardeenlll970l ; 
lThorne|[l97i ). iKing et al. I (|2005h argue instead that most 
MBHs have very low or no spin, due to preferential accre- 
tion of co unter-rotating m aterial, or to short-lived accretion 
episodes (|King et al.ll2005h . Waiting for additional observa- 
tion^, we consider here two extreme cases, that likely allow 
us to bracket the typical configurations: either that all MBHs 
have exactly null spin, or that all MBHs have a = 0.9. The 
latter value is slightly lower than the canonical a — 0.998 
jThornd 11971 ). but it is consistent with magnetohydrody - 
namical simulations of disc accretion (|Gammie et al.ll2004T ). 
We also run a control simulation where we set the recoil 
velocity to zero for all MBH mergers. 

For every galaxy merger we track jointly the dynamical 
evolution of the MBHs and of the host halo. In addition to 
the dynamics of MBH binaries, as described above, we trace 
the mass evolution of their hosts, including the increase due 
to galaxy mergers and the decrease due to mass stripping of 
the halo within the gravitational potential of the halo. Our 
treatment is very simple: we integrate the equation of motion 
of the satellite in the gravitational potential of the cluster 
(assuming a non singular isothermal sphere), including the 
dynamical friction term. At every step of the integration we 
compare the density of the satellite to the density of the 
cluster halo at the location of the satellite. Tidal stripping 
occurs at the radius within which the mean density of the 
satellite exc eeds the density of the galaxy interior to its or- 
bital radius (|Tavlor fc BabufeOOll ). We trace tidal stripping 
of all satellites from the time of the merger to z = 0. 



3 BLACK HOLES OCCUPATION FRACTION 
AND MASS SCALING 

By applying the assumptions described in Section 2 to the 
assembly history of galaxies in a region that can be iden- 



3 The spin is a measurable parameter, wi th a very high accuracy , 
in the gravitational waves LISA signa l llBarack fc Cutler] I200I 
iBerti et alj|2005l ; lLang fc Hughes! |2006| ; IVecchioll2004) . Gravita- 
tional waves emitted when a stellar mass MBH, or even a white 
dwarf, or neutron star, falls into a massive black hole can be used 
to map the spacetime of MBHs, and as a consequence, MBH spins. 
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Figure 1. Occupation fraction of MBHs in Virgo galaxies. 
Hatched histogram: Decarli et al. 2007. Dashed histogram: no 
recoil. Dotted histogram: a = 0. Solid histogram: a = 0.9. From 
bottom- left panel, clockwise: model I (higheff), model II (mid- 
eff), model III (peak3), model IV (peak3.5). Model histograms 
are based on total baryonic masses. 



tified as a cluster at z = 0, we can determine the fraction 
of the galaxies which host a MBH, and if there are trends 
with the size and mass of the host galaxy. We note here 
that when we refer to the mass of the host, we assume for 
simplicity the total baryonic mass, which we scale with the 
cosmic baryon fraction (14%). The stellar mass M* of galax- 
ies, as well as the dynamical mass within the optical radius, 
Mgai, are undoubtedly non constant functions of the halo 
mass, environmental density, and star formation history of 
the satellite itself. Total baryonic masses are therefore up- 
per limits to the measurable stellar and dynamical masses 
of the hosts. Given our assumptions on MBH formation and 
growth, the occupation fraction does not sensibly depend on 
bulge mass, but rather on the total mass of a galaxy, consid- 
ered a tracer of the cosmological formation time (linked to 
the halo bias), and of the dynamical evolution (within the 
cluster potential). We stress, however, that pure disc galax- 
ies are likely unsuitable for MBH formation (see discussion 
in section 2). 

We determine the BHOF as a function of halo mass 
(compare with the sta ndard definition for galaxies, e.g., 
iPeacock fc Smithll2000D . which is shown in Figure [T]for the 
four different models I, II, III, IV. We compare our BHOF to 
the di stribution of AGN in Virgo galaxies bv lDecarli et al.l 
<|2007h . who find signs of AGN activity in all late type 
galaxies in Virgo with galaxy mass above 3 x 10 11 Mq. 
Clearly the AGN occupation function is a lower limit to the 
BHOF, that is at every galaxy mass: BHOF(M 9ai ,AGN)^ 
BHOF(M 9a ,Ki- 

The BHOF of our models where the recoil velocity is 
set to zero informs us of the "nature" scenario, in that only 
the initial conditions of MBH formation at high redshift en- 
ter into the final BHOF that we find at z = 0. We see that 
in all models I, II, III and IV, the BHOF is of order unity 



Figure 2. MBH mass against baryonic mass. Circles: no recoil. 
Triangles: a = 0. Stars: a = 0.9. From bottom-left panel, clock- 
wise: model I (higheff), model II (mideff), model III (peak3), 
model IV (peak3.5). 



only above a certain galaxy mass. It is apparent that the 
existence of a threshold in galaxy mass (as traced also by 
the velocity dispersion) is a natural outcome of MBH forma- 
tion scenarios where the formation redshifts are pushed to 
high redshift. We note here that it is indeed necessary that 
MBHs are formed early on, as luminous quasars are detected 
at redshift 5 and higher (e.g., Fan et al. 2004). Obviously, 
the lower the typical MBH formation redshift, the higher 
the fraction of dwarf galaxies that can host a MBH. Even 
without additional physical processes, our analysis shows 
that the possible decrease in BHOF at low galaxy mass is 
an expected feature in models that are based on realistic 
MBH seed formation mechanisms. A similar behavior, that 
is, an increasing BHOF with galaxy mass, is expected also 
in field galaxies (I Volonteri et al.ll2007t ). but it is exacerbated 
in cluster satellite galaxies, whe re coalescence rounds (using 
the terminology of ISesanal I2007D are truncated. 

The dynamical evolution of the MBH population mod- 
ifies the BHOF in the direction of decreasing it, as MBHs 
are prone to ejections, and the ejection probability depends 
on the galaxy mass, via the escape velocity. Although the 
changes are not dramatic, that is MBH k icks do not de- 
plete most z = galaxies of their MBHs (|Volonteril 120071 : 
ISchnittmanll2007l '). the effects are non negligible, especially 
for galaxies with mass below 10 11 M . The BHOF becomes 
a steeper function of the galaxy mass, and the threshold 
for BHOF close to unity shifts toward s higher masses. If w e 
compare our models to the results bv lDecarli et al.l (|2007t ) . 
who find that all galaxies with mass above 3 x 10 11 Mq show 
signs of AGN activity, hence host a central MBH, we can rule 
out very high biased models (model IV, peak3.5) with high 
recoil velocity (solid lines in Figure [1]), so if MBH forma- 
tion is indeed very biased, either MBH spi ns align during the 
pre-m erger orbital decay, as suggested by iBogdanovic et alj 
(120071 1. or MBHs have small spins. 

An additional prediction that arises from our models is 
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that galaxies which have experienced sustained tidal strip- 
ping in the deep potential of the cluster host MBHs which 
lie above the expected Mbh — M ga i correlation. After a 
merger episode that grew MBH and galaxy, if the orbital 
decay of the host brings it in dense regions, tidal stripping 
will decrease the galaxy mass, while the MBH mass will 
not be modified. The effect is more pronounced if MBHs 
are widespread, as shown in Figure [2] If MBHs form at 
early times, and evolve jointly with their hosts, with accre- 
tion bringing the MBH s exactly on the Mbh — a* relation 
(|Di Matteo et al.1 12005] ). then this population of "overmas- 
sive" MBHs is unavoidable in galaxy clusters. We strongly 
stress that our results are robust in a qualitative way, al- 
though they should not be assumed quantitatively correct. 
For instance, the slope of our Mbh — M ga i correlation is 
of order 5/3, because of our assumed rigid scaling between 
halo properties and MBH properties. Regardless of the ex- 
act slope, however, the presence of stripped galaxies with 
overmassive black holes is a robust prediction of models 
where MBHs form, and grow, early on. We also draw the 
attention to the possible environmental dependence of the 
Mbh — Mgai correlation, but not of the Mbh — cr*, as the 
central velocity dispersion is not likely to be affected by the 
stripping. 

There is observational evidence of overm a ssive MBHs 
for a given M* in Virgo. iKormendv et all (1997) used 
ground-based spectroscopic observations of NGC 4486B, the 
satellite galaxy to NGC 4486 (M87) , along with isotropic dy- 
namical models to find a central dark object of mass M — 
6^2 x 10 s Mq , though three-integr al, axisymmetric models 
cannot rule out absence of a MBH. IMagorrian et al l (| 19981 ) 
later modeled NGC 4486B with two-integral, axisymmetric 
dynamical models to find a mass oi M = 9.2+Q Q33 x 10 9 M . 
The galaxy contains a double nucleus, which may pose prob- 
lems for isotropic and axisymmetric models, but the pres- 
ence of a double nucleus, itself, is explained by a n eccentric 
disk of stars, which requires a MBH as in M31 (|Tremainel 
Il995l ) So while the mass estimate may be uncertain, there 
is strong evidence for a MBH in NGC 4486B. Even if the 
mass estimate is high b y a factor of a few, wi th a stellar 
mass of ~ 6 x 10 9 M (|Kormendv et alj 1 19971) the MBH 
is more massive than those in galaxies wit h comparable 
bulge mass (see, e.g.. IMagorrian et al.lll998l ). NGC 4486B 
is almost certainly tidally stripped, though this could be 
primarily a result of mu ch stronger interactions with M87 
((Evstigneeva et al]|2007l ) rather than milder stripping from 
the cluster potential. Another potential example in Virgo is 
the SO galaxy NGC 4342 with a measured black hole mass of 
M = 3.0±] l x 10 s M B and bulge mass of M = 1.2 x 10 10 M Q 
|Cretton fc van den Bosch] 1 19991 ). m aking it another ex - 
treme high outlier by a factor of ~ 30 (|Haring fc Rixll2004l ). 

We note that the population of "undermassive" MBHs, 
signature of th e gravitational r ecoil, described, e.g. in 
IVolonteril (|2007l ); IVolonteri et ail l|2007l ) is basically absent 
here, due to the more simplistic modeling we performed. We 
have in fact assumed that whenever two galaxies with MBHs 
merge, with a mass ratio larger than 1 : 10, the MBHs are 
placed directly on the Mbh — o* ■ This consequently puts 
MBHs preferentially on the Mbh — M ga i correlation. The 
cases where the MBH is below the mass predicted by the 
Mbh — Mg a i are those in which the galaxy has not had any 
major merger recently, so the MBH mass has not grown, 



while the galaxy mass has (via, minor mergers, cfr. the dis- 
cussion m IVolonteri et~ai1l2003l ). 



4 CONCLUSIONS 

We explored the influence of formation epoch of MBHs, bias 
of their hosts, and MBHs dynamical evolution on the occu- 
pation distribution of MBHs in galaxy clusters. The possi- 
ble decrease of the occupation f ractio n at l ow galaxy mass 
propo sed by IWehner fc Harris] (|2006l ) and IFerrarese et alj 
(2006) is not a surprising result and follows naturally from 
the evolution of the MBH population of galaxies in clusters. 

• When the formation mechanisms of MBH seeds are 
taken into consideration, implying formation of MBHs in 
massive high redshift halos, the BHOF in cluster galaxies in 
an increasing function of galaxy mass, in line with observa- 
tions of the AGN fraction in Virgo. 

• The exact mass threshold above which BHOF=l de- 
pends on the details of the formation mechanism (host 
masses and redshift of formation, "nature") and on the dy- 
namical evolution, including MBH spin magnitude ("nur- 
ture"). 

• The repercussions of "nurture" are magnified in cluster 
galaxies. If a satellite galaxy looses its MBH due to a dy- 
namical interaction, it has a negligible chance of capturing 
a new one following a subsequent galaxy merger, except for 
the central galaxy in the cluster. 

• We also predict that if MBHs co-evolve with galax- 
ies during galaxy mergers, and satellite galaxies experience 
tidal stripping during their orbital evolution, then a popu- 
lation of galaxies where the MBH mass lies above the stan- 
dard Mbh — Mgai is expected. Two such examples could be 
NGC 4486B and NGC 4342. 

We emphasize that our models do not in any way imply 
that MBHs and nuclear clusters are mutually exclusive. We 
however predict that the occupation fraction decreases with 
galaxy mass. If, as suggested by Wehner & Harris (2006), 
lacking a MBH is a necessary condition for nuclear cluster 
formation, radio or hard X-ray observations will inform us of 
the complementary or mutually exclus ive essence of MBHs 
and nuclear clusters. iGallo et alj l|2007l ) recently observed 32 
galaxies in Virgo with the Chandra X-r a y Obs ervatory, and 
found, in agreement with iDecarli et all ((20071 ) that nuclear 
X-ray activity increases with the mass of the host galaxy. 
Intriguingly, at least in one case, VCC1178, a nuclear X- 
ray source is det ected jointly with a central star cluster. 
ISeth et alj j|2006l ) also report that a small AGN might be 
hosted within the core of the nuclear star cluster in NGC 
4206. 
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